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Motivated by the observation of the decay B — > K*"f by CLEO, we have systematically 



Abstract 



analyzed the two-body weak radiative decays of bottom and charmed hadrons. There exist 



two types of weak radiative decays: One proceeds through the short- distance b — > S7 transi- 
tion and the other occurs through PU-exchange accompanied by a photon emission. Effective 
Lagrangians are derived for the ^-exchange bremsstrahlung processes at the quark level and 
then applied to various weak electromagnetic decays of heavy hadrons. Predictions for the 
branching ratios of B° — ► -D*°7, A° — > £^7, 5° — * S°7 and 2° — > S' c 7 are given. In par- 
ticular, we found B(B° — > D* 0, ~f) « 0.9 x 10~ 6 . Order of magnitude estimates for the weak 
radiative decays of charmed hadrons: D° — > K*°j, A+ — > £+7 and H° — > H°7 are also 
presented. Within this approach, the decay asymmetry for antitriplet to antitriplet heavy 
baryon weak radiative transitions is uniquely predicted by heavy quark symmetry. The elec- 
tromagnetic penguin contribution to A° — > A7 is estimated by two different methods and its 
branching ratio is found to be of order 1 x 10~ 5 . We conclude that weak radiative decays of 
bottom hadrons are dominated by the short- distance b — > 57 mechanism. 



I. Introduction 

The recent observation of the weak radiative decay B — ► K*^y by CLEO [1] with the 
branching ratio (4.5 ± 1.5 ± 0.9) x 1CT 5 confirms the standard-model expectation that this 
decay mode is dominated by the short-distance electromagnetic penguin transition b — ► s'y. 
Naively, it is tempting to think that B — ► D*'-/ will be the dominant weak radiative decay 
of the B meson as it is not suppressed by quark mixing angles. However, owing to the 
large top quark mass, the amplitude of b —* s'-f is neither quark mixing nor loop suppressed. 
Moreover, it is largely enhanced by QCD corrections. As a consequence, the short-distance 
contribution due to the electromagnetic penguin diagram dominates over the W -exchange 
bremsstrahlung. This phenomenon is quite unique to the bottom hadrons which contain a 
heavy b quark; such a magic short-distance enhancement does not occur in the systems of 
charmed and strange hadrons. For example, it is known that the mechanism s — > d'y plays 
only a minor role in the radiative decays of kaons and hyperons. 

In Ref . [2] we have systematically studied the flavor-conserving electromagnetic decays of 
heavy mesons and heavy baryons. Various photon coupling constants are related through the 
usage of heavy quark symmetry. For example, the B*B* r y coupling, which is very difficult 
to measure in practice, is related to the B*B , ~f coupling via heavy-quark spin symmetry. 
The coupling constants appearing in the Lagrangians depend only on the light quarks and 
can be calculated in the nonrelativistic quark model. Consequently, the dynamics of the 
electromagnetic transitions for emissions of soft photons and pions is completely determined 
by heavy quark and chiral symmetry, supplemented by the quark model. The purpose of the 
present paper is to extend our previous work to the weak radiative decays of heavy hadrons. 

At the quark level, there are three different types of processes which can contribute to the 
weak radiative decays of heavy hadrons, namely, single-, two- and three-quark transitions 
[3]. The single-quark transition mechanism comes from the so-called electromagnetic penguin 
diagram. Since the penguin process c — » wy is very suppressed, it plays no role in charmed 
hadron radiative decays. We will thus focus on the two-body radiative decays of bottom 
hadrons proceeding through the electromagnetic penguin mechanism b — > 57: 

A° -> S° 7 , A° 7 , ~° -> S° 7 , S b " -> S" 7 , fi 6 " - ^7- l ' ' 

There are two contributions from the two-quark transitions: one from the VT-exchange dia- 
gram accompanied by a photon emission from the external quark (see, for example, Fig. 1), 
and the other from the same VF-exchange diagram but with a photon radiated from the W 
boson. The latter is typically suppressed by a factor of m q k/M^r {k being the photon energy) 
as compared to the former bremsstrahlung process [4]. For bottom hadrons, the dominant 



decays which occur through the quark-quark bremsstrahlung bd — > cwy or bu — ► cd'-f are: 

(1 21 

where we have followed the convention that a 5 meson contains a 6 quark and that Sq (Hq) 
denote antitriplet (sextet) heavy baryons. For charmed hadrons, the Cabibbo-allowed decay 
modes via cu — »■ SJ7 or erf — > MS7 are 

A+ -> S+7, ~° - ~° 7 . l • J 

Note that some decay modes in (1.1) also receive contributions from VF-exchange bremsstrahlung, 
but they are suppressed by quark mixing angles. Finally, the three-quark transition involv- 
ing VF-exchange between two quarks and a photon emission by the third quark is quite 
suppressed because of very small probability of finding three quarks in adequate kinematic 
matching with the baryons [3,5]. 

To summarize, the two important mechanisms for weak radiative decays of heavy hadrons 
are PF-exchange bremsstrahlung and the electromagnectic penguin transition b — ► s'y. Since 
the effective Lagrangian for the latter is known, the calculation for the radiative amplitude 
induced by the penguin diagram appears easier at first sight. 

The W -exchange bremsstrahlung effect is usually evaluated under the pole assump- 
tion; that is, its amplitude is saturated by one-particle intermediate states. When dealing 
with weak radiative decays of heavy hadrons, one encounters two predicaments. First, the 
hadronic matrix elements for the processes (1.1) are evaluated at q 2 = for a real photon 
emission, whereas heavy quark symmetry and the quark model are known to be more reli- 
able at zero recoil kinematic point where q 2 is maximum. (The quark-model wave functions 
best resemble the hadron states in the frame where both hadrons are static.) Second, the 
intermediate states appearing in the pole diagrams for the processes (1.2) or (1.3) are very 
far from their mass shell. For example, the four-momentum squared of the D pole in the 
decay B — > _D*7 is m 2 B (see Fig. 3). This means that the residual momentum of the D meson 
defined by P M = m D t),, + k^ must be of order rriB, so the approximation k/m D << 1 required 
by the heavy quark effective theory is no longer valid. Also, the quark-model prediction for 
the photon coupling constants is presumably reliable only when both hadrons are nearly on 
their mass shell. The question is then how to extrapolate the hadronic matrix elements from 
zero recoil to maximal recoil, and the photon couplings from the on-shell point to off-shell? 
In principle, one can treat the intermediate state as an on-shell particle and then assume that 
off-shell effects of the pole can be parametrized in terms of form factors. Such form factors 
are basically unknown, though they are expected to become smaller as the intermediate state 



is more away from its mass shell due to less overlap of initial and final hadron wave func- 
tions. Consequently, based on heavy quark symmetry and the nonrelativistic quark model, 
at best we can only predict the upper bound of the decay rates for the radiative decays in 
the category of (1.2). 

We will present in this paper a different but more powerful approach to the W^-exchange 
bremsstrahlung processes. The fact that the intermediate quark state in these processes is 
sufficiently off-shell (see, e.g. Fig. 1) and the emitted photon is hard suggests the possibility 
of analyzing these processes by perturbative QCD. As a first step in this direction, we study 
the tree amplitudes responsible for these processes and derive a gauge invariant effective 
five-point interaction for the quark-quark bremsstrahlung bu — > cd'-f or bd — > cu'j. 

The physical mass of a heavy hadron differs from the heavy quark mass by an amount 
of order Aqcd- This difference is due to the presence of the light quark(s). It is therefore 
reasonable to assign a constitutent mass of order Aqcd to the light quark(s) inside a heavy 
hadron. In addition, the light quarks move, on average, with the same velocity as the heavy 
quark. We will make the simplifying assumption of neglecting the relative Fermi motion. 
Thus, the heavy quark and the light quark(s) in a heavy hadron move with equal four velocity. 
This momentum parametrization has the advantage that the resulting effective interaction 
is local and manifestly gauge invariant. In Sec. Ill we will show explicitly for the meson case 
that the effective Lagrangian and the pole model approaches are indeed equivalent, but the 
former is much simpler and provides information on the form factors. 

Armed with the effective Lagrangian for the W^-exchange bremsstrahlung, we are able to 
study various radiative decay modes of bottom and charmed hadrons listed in (1.2) and (1.3), 
bearing in mind that this approach presumably works better when both the initial and final 
hadrons contain a heavy quark. We will use the factorization method, which is known to work 
well for nonleptonic weak decays of heavy mesons, to evaluate the mesonic matrix elements. 
As for the baryon radiative decays, we will demonstrate that heavy quark symmetry leads 
to a nontrivial prediction for antitriplet to antitriplet heavy baryon transitions: The ratio 
of parity-conserving and parity- violating amplitudes is uniquely predicted. Baryonic matrix 
elements will be calculated using the MIT bag model. 

This paper is organized as follows. Local effective Lagrangians for the quark-quark 
bremsstrahlung processes are derived in Sec. II. We then apply this approach to weak 
radiative decays of heavy mesons in Sec. Ill and to various bottom decays in Sec. IV. 
Discussion and conclusion are presented in Section V. Some preliminary results have been 
reported earlier [6]. 



II. Effective Lagrangians for Weak Radiative Decays 



In this Section we will present the effective Lagrangians for the penguin transition b — ► 
s+7 and for the W^-exchange bremsstrahlung processes b+u — > c+d+7 and b+d — > c+w+7. 



The effective Lagrangian for the short-distance 6 
tions reads [7] 
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where a ■ F = a llu F^ v , and we have neglected contributions which vanish for a real photon 
emission, V^- is the quark mixing matrix element, F 2 = F 2 (x t ) — F 2 (x c ) = F 2 (x t ) with 
Xi = m%/Myy, and 
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It is easily seen that F 2 is a smooth function of the top quark mass. For m t = 150 GeV 
and Aq CD = 200 MeV, we find F 2 (x t ) = 0.34 and F 2 {x t ) = 0.73, so that the radiative decay 
b — > s'-f is enhanced by QCD corrections by a factor of 2. The radiative decays B — »■ K*^ 
and I? — ► 07, mediated by this penguin mechanism, have been studied extensively in the 
literature. In Section IV we will apply the effective Lagrangian (2.1) to the decay A5 — > A7. 

We next turn to the W-exchange bremsstrahlung processes bd — > cu'-f and bu — > cd'-f. The 
difficulty associated with highly off-shell intermediate states mentioned in the Introduction is 
easily overcome at the quark level. The propagator of the highly virtual quark can be reduced 
to a constant by energy-momentum conservation. The above photon emission reactions are 
then described by an effective five-point local interaction which is also gauge invariant. To 
begin with, we note that the relevant QCD-corrected effective weak Hamiltonian is given by 
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where ((71(72) = 9i7a»(1 — 75)92- The Wilson coefficient functions c±, evaluated at the scale 
/i = mb, have the values 

c+(m 6 ) ~ 0.85 , c-{m b ) ~ 1.38 . (2.7) 

We first consider the photon emission process &d — ► cu'j. The amplitudes mediated by 
the operator Oa (see Fig. 1) are 

At = ee c u c Y , u 7^(1 - lb)u b Vdlu{^ - 7s)f u , (2.8a) 

Pc ~T If m c 

A 2 = ee b u c Y(l - 75)7 u ^UbVd^il - 7 5 )v u , (2.86) 

Pb lp m b 

A 3 = ee d u c j u (l - 7 5 )m^7 m — - — —. 7„(1 - 75)^, (2.8c) 

— FJ + r ~ m d 

A 4 = ee tt w c 7"(l - 75)«6^7^( 1 ~ 7s)— 7 « 7^«, (2.8d) 

-?« - f - m„ 

where A; is the photon momentum. 

We will parametrize the quark momenta in terms of velocities; this is more suitable when 
dealing with heavy quark symmetry: 

p b = m b v, pa = m d v, p c = m c v', p n = m u v'. (2.9) 

Since the b quark is heavy, we have set the velocity of the d quark to be the same as the b 
quark so that they will move together to form a bound meson state. Likewise, v a = v c = v '. 
At this point we wish to emphasize that the light quark masses appearing in (2.8)-(2.9) are 
of the constituent type. This is attributed to the fact the typical Fermi momentum of the 
quarks in a hadron is of order Aqcd- Consequently, although the current quark masses of 
the light u and d quarks are only of order 10 MeV, their off-shellness is of order Aqcd- We 
thus choose to have the light quarks close to their mass shell, so that p q ~ m q v with v 2 = 1 
and m q being the constituent quark mass. Obviously, this parametrization (2.9) does not 
provide a complete description of the Fermi motion inside the bound state. Nevertheless, 
it does take into account its average effect by giving a constitutent mass of order Aqcd to 
the light quarks. This parametrization greatly simplifies the calculation by eliminating the 
photon's coupling to the convection currents and making the effective interaction local and 
manifestly gauge invariant. 

With the momentum parametrization given by (2.9), it is easily seen that the contribu- 
tions from the convection current add up to zero and the amplitude arises entirely from the 



magnetic moments of the quarks: 
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where rrii = m b + m d and rrif = m c + m u . The amplitude corresponding to Fig.2 induced 
by the operator Ob can be obtained from Eq.(2.10) by the substitution u b <-► v u . The above 
amplitudes can be further simplified by considering the commutator and anticommutator 
relations 

2{<V» 7a} = e^ A a7 Q 75, (2.11a) 
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We find that for the photon emission process bd — > cu'y, we can simply replace the operator 
0± in (2.4) by 0± so that the effective Hamiltonian is given by 

W*(M -> czZ 7 ) = ^jf cd V: d {c + F + + c_OF), (2.12) 
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Of = c 7 "(l - 7 5 )MV(1 - 7s)w ± c 7 M (l - l 5 )udY(l - -y 5 )b, (2.14a) 
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Similarly, for the ^-exchange bremsstrahlung bu — > 007, we have 

0£(&u-cd 7 )=^ { (e c ^-e^)(^ + ^) 
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where now rrii = m b + m u , to/ = m c + m d . 



(2.15) 



Strictly speaking, terms in (2.13) and (2.15) proportional to 1/to c and l/m b should have 
been dropped since we have not included corrections of the same order to those proportional 
to l/m u or l/m d . These l/ro c and l/m b terms are kept mainly for phenomenological reasons. 
We encounter a similar situation in the electromagnetic decays of heavy hadrons [2,8] in which 
the common practice is to retain the 1/toq contributions due to the heavy quark's magnetic 



moments but to neglect the I/tuq corrections to the matrix elements of the light quark's 
electromagnetic currents. Numerically, the l/m c terms are especially significant because m c 
is not very large and the c quark carries 2/3 units of charge. The effective Lagrangians (2.13) 
and (2.15) are the main results in this section. We will apply them and (2.1) to the weak 
radiative decays of heavy mesons and heavy baryons in Sections III and IV, respectively. 

III. Applications to Heavy Meson Decays 

As shown in the Introduction, the radiative decay modes of interest for B mesons are 
B —>■ K*j, B s — > 07 which receive short-distance contributions from the electromagnetic 
penguin b — > s'j transition, and B —>■ D*^ proceeding through the W exchange accompanied 
by a photon emission. The general amplitude of weak radiative decay with one real photon 
emission is given by 

A[B(p) -* P*(q)i(k)] = ie^eWe^qPfiik?) 

+ ^[e;(m|-m^)-(p + g)^-fc]/ 2 (fc 2 ), 

where e and e* are the polarization vectors of the photon and the vector meson P*, respec- 
tively, the first (second) term on the r.h.s. is parity conserving (violating), and k 2 = 0. The 
decay width implied by the amplitude (3.1) is 

r(B -> P* T ) = -^ ^ ~f P * ? (l/il 2 + 4|/ 2 | 2 ). (3.2) 
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3.1 Effective Lagrangian approach for B" — > D* u j 

Since the radiative decay B — > K*j has been discussed extensively in the literature, we 
will only focus on the second-type mode, namely B° — > D* 0/ ~f. Our goal is to see if the 
tree-level ^-exchange with a photon emission is comparable with the short-distance b —>■ s'j 
mechanism. We shall use the factorization method (for a review, see Ref.[9]) to evaluate the 
hadronic matrix elements. It follows from Eq.(2.12) that 

A(B° -► D*%) = -(D*°-f\H cfi \B ) 

-%V cb V: d a 2 -^(D* \cj,(l - l5 )u\0)(0\d lu (l - 75)&|S°) 
v i / 
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«2 = ^{c- -c+), (3.4) 

^n the conventional vacuum insertion method a2 is equal to (2c+ — c_)/3, while it is (c+ — c_)/2 in the 



and m>i = mb + ma ~ tub, rrif = m c + m u « mo*- The one-body matrix elements appearing 
in (3.3) have the expressions 

(0\Ap\P(p)) = ifp Plt , (3.5a) 

(0\Vp\P*(p,e*))= if v mp*s;. (3.56) 

Therefore, 
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With the substitution of (3.6), (3.3) has the desirable amplitude structure indicated by (3.1). 
We find 
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Substituting (3.7) into (3.2) gives the decay rateQ 



r(B°-Z>*° 7 ) = |G^,|WlS.m B 1 



c_ - c A 



X 



(e^ (^ + X\ + e 2 d ( ^ + -\) + 2e u e d ^ (-^- + -± 



In order to have a numerical estimate, we adopt the following mass parameters 

m B o = 5279.0 MeV, m D *o = 2006.7 MeV, 
m u = 338 MeV, m d = 322 MeV, m s = 510 MeV, 



(3.9) 



(3.10) 



from the Particle Data Group [15], and m c = 1.6 GeV, m& = 5.0 GeV. Using /d* = 200 MeV, 
f B = 190 MeV, V cb = 0.040 [16], and t(B°) = 1.50 x 10~ 12 s [15], we obtain the branching 
ratio 

B(B° -> £>*%) = 0.92 x 10~ 6 . (3.11) 

large N c approximation in which the Fierz-transformed contributions characterized by the color factor 1/N C 
are dropped [9]. The leading 1/N C expansion is known to work well for nonleptonic weak decays of charmed 
mesons. In bottom meson decays, the magnitude of a-i determined from the measured B — > ipK, tpK* 
rates is in agreement with that predicted by the large N c approach [10-12]. However, contrary to what 
expected from the same approach, the sign of a-i is found to be positive by recent CLEO measurements of 
B — > D-k, Dp, D*ir, D*p decays [11-13]. Thus we take a 2 to be that given by (3.4). 

2 The result (3.9) was also obtained by Mendel and Stiarski [14] in a different approach except for the 
Wilson factor ai being replaced by (2c + — c_)/3 in the latter work. 

3 The values of the constituent quark masses are given on p. 1729 of Ref.[15]. 
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It is evident that the weak radiative decays of the B mesons is dominated by the electro- 
magnetic penguin diagram. The suppression of B — > D*^ relative to B — > K*^ is mainly 
attributed to the smallness of the decay constants Jd* and fp occurred in weak transitions. 

In the factorization approximation the three decay amplitudes B° — > D* 0/ -f, B*° — > D* 0, y 
and B*° — »■ D 0/ ~f are all related since fp = fy by heavy quark symmetry. However, the 
branching ratios for the latter two reactions are expected to be much smaller than that for 
B° — ► D* 0/ ~f as B*° has a dominant electromagnetic decay. 

3.2 Pole model approach for B° — ► Z}*°7 

Before proceeding further, we would like to compare our present formalism with the 
conventional long-distance pole mechanism which has been applied to the decay B — ► K*^ 
before [17]. Note that the intermediate B state is absent in Fig. 3d as the BB^ coupling is 
prohibited. If we just focus on the low-lying intermdiate states in the pole diagrams depicted 
in Fig. 3, the amplitudes of the first three pole diagrams are 

M a = (D*°(q)\C cm \D°(p)} 2 l 2 (D (p)\C w \B°(p)}, (3.12a) 

m B — m D 

M b = (D*°(q)\C w \B*°(q)) 2 l 2 (B*°(q)\£ cm \B (p)), (3.126) 

m D * — m B * 

M c = (D*°(q)\£ em \D* (p)) 2 1 2 (D*°(p)\C w \B°(p)), (3.12c) 

m B — m D * 

where the amplitudes M a and Mb are parity conserving, while M c is parity violating. Since 

the intermediate pole states are far from their mass shell, we write the photon couplings in 

Fig. 3 as 

g DD *-M 2 = m B) = 9 D (q 2 )g DD ,^ 

9 b *bM 2 = w2 d) = 9 B *(q 2 )g B * Bl , (3-13) 

9 D * D «M 2 = m2 B) = 9 D M 2 )g D * D *-,i 

where g DD * , 9 B , Bl and g D * D * are on-shell photon coupling constants which can be calculated 
in the nonrelativistic quark model. In Eq.(3.13), g D , g B , and g D , are form factors accounting 
for off-shell effects. They are normalized to unity when mesons are on shell; for example, 
9 D ijn 2 D ) = 1, g B *{m 2 B ,) = 1. 

For P,P* = (Qq), we find from Eqs.(2.19) and (2.30) of Ref.[2] that 

(D*°(q)\£ em \D°(p)) = ^Kls^e^fc"^ ', (3.14a) 

(B*°(q)\C cm \B°(p)} = ig B ^m 2 D )g Bt0B0i e, va ^k v p a e^, (3.146) 

(D*°(q)\jC em \D*°(p)) = g D *(m 2 B )g D , 0D ^(k»e u - fc^(?)e w (pK.. (3.14c) 

10 



Because the intermediate state is far from its mass shell, we have introduced the form factors 
defined in Eq.(3.13). In the heavy quark effective theory, the on-shell P*P'~f and P*P*7 
coupling constants are related to each other and are given by [2] 

g PP .~ = -2J—^-(e q d + e Q d') = ej—^(e q (3 + e Q (3'), (3.15a) 



9 P * Pl = -2J—^(e q d + e Q d') = eJ—^(e q P + e Q [3'), (3.156) 

y mp* y mp* 

9p. P .^ = ~2(e q d - e Q d') = e(e q (3 - e Q (3'), (3.15c) 

where e q is the charge of the light quark q (not q), and cq is the charge of the heavy quark 
Q. The coupling f3' (or d') is fixed by heavy quark symmetry to be 1/rriQ, while j3 (or d) is 
independent of heavy quarks and cannot be determined by heavy quark symmetry alone. In 
the constituent quark model, /3 is given by [2] 



(5= —. (3.16) 

m„ 



Therefore, 



e, 



W,-(^ + ^), (3,76) 

^^ = e (S:"^)' (3 - 17c) 

where the small difference between mp and mp* has been neglected. Note that it is important 
to keep the contribution from the magnetic moment of the charmed quark since it is not 
particularly heavy. In general, it is expected that the form factors appearing in Eq.(3.14) 
become smaller as the hadron is more away from its mass shell owing to less overlap between 
initial and final meson wave functions. 

Using the factorization method, we obtain 

(D°(p)\jC w \B°(p)) = ^L Vcb V: d (c- - c + )f D f B m%, (3.18) 

and use has been made of p 2 = m\. Likewise, we obtain 

(D*°(q)\£ w \B*°(q)) = ^V^V^c. - c + )f D f B m D m B (e* B ■ e* D ), (3.19a) 

(D*°(p)\C w \B°(p)} = -^= Vcb V^(c.-c + )f D f B m D *(p-s* D ). (3.1%) 
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Putting everything together and using the relation 

A 

we finally find 



m z 



J D* 



Til r> / 9 \ 



m D m B 



-9^onoJ B *(m 2 D ) 



"B "°D 



" l B* " l D* 



fo 



/ 2 (0) = -^9 fl ,c D ,o 7 9 D .K). 



Comparing (3.21) with (3.7) yields 



# D ( m l) 



m 



_u 



m 



m 



2 • 



B . {™DJ 



9 D A m B) 



D* 



(3.20) 

(3.21a) 
(3.216) 

(3.22) 

i B IIL B - ,li D , 

where use of the heavy quark symmetry relations m B * = tub, vri£,* = mo and fp* = fp 
has been made. We note that in the effective Lagrangian approach there is an additional 
term proportional to (e^/m^ — eb/mb) in f'2. What is the counterpart of this term in the 
pole model? Evidently, it must come from a p-wave 1 + B resonance state (see Fig. 3d): 
The El B° (1 + ) B° (0~)7 transition coupling is proportional to {e^/ raj — eb/mb) provided that 
B(l + ) is ap-wave spin-singlet, while -B(l + ) — D* weak transition is parity violating. We thus 
see that both approaches are consistent with each other. However, the effective Lagrangian 
approach is simpler and it also provides the information on the form factors, as shown in 
(3.22). 



m, 



in 



2 ' 



IV. Applications to Heavy Baryon Decays 

In this section we will first focus our attention on the short-distance penguin effect in the 
decays At — > £7 and A& — > A7 and then turn to the weak radiative decays of the antitriplet 
bottom baryons, namely S° — »■ E^(S' C )7 and A° — > S°7. Since the weak radiative decay of 
B mesons is dominated by the short-distance b — > 57 transition, it is natural to expect that 
the same mechanism also works for bottom baryons. The general amplitude of baryon weak 
radiative decay reads 

A(Bi -► B f i) = iu f (a + ln 5 )a^k v u h (4.1) 

where a and b are parity-conserving and -violating amplitudes, respectively. The correspond- 
ing decay rate is 

T(B t - B n ) = ±- (Hi^l) 3 ( | Q | 2 + \bf). (4.2) 

8vr \ mi J 

4.1 Penguin induced baryon radiative decays 
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The electromagnetic penguin transition 6 — »■ s + 7 at the quark level is well understood in 
QCD, and its effective Lagrangian is given by (2.1). However, its hadronic matrix elements 
depend on the nature of the bound states, and there are no known methods to calculate them 
from first principles. In this Section, we will use two different methods to estimate the rates 
for the decay A5 — *• A + 7. In the first method, both the 6 and s quarks are treated as heavy. 
Heavy quark symmetry then gives model independent prediction. Although we are able to 
compute the correction of order m s /mb and l/m s , higher order corrections are needed for 
a realistic comparsion with experiment. They are in principle well defined in QCD; but it 
is impossible to compute them all. In the second method, only the b quark is treated as 
heavy so the form factors of tensor currents needed are related by heavy quark symmetry to 
those of vector and axial vector currents, which are then evaluated in the MIT bag model. 
The MIT bag model includes the m s effects to all orders, but the reliability of the model is 
hard to assess. It is clear that neither of the two approaches is very satisfactory. Together, 
however, they provide a more or less consistent order of magnitude estimate of these weak 
radiative decays. 

In the first method we treat the s quark as a heavy quark and then take into account the 
l/m s and QCD corrections. Despite that the effective mass of the s quark in hyperons is 
only of order 500 MeV, it is not small compared to the QCD scale and we thus expect to see 
some vestiges of heavy quark symmetry. In the heavy s quark limit, the hyperon A behaves 
as an antitriplet heavy baryon B%, while S° as a sextet baryon Bq. From Eq.(2.1) we obtain 



A(A° b -► hyperon + 7 ) SD = %% -£?F 2 V tb V? s m h e»k» 

x <hyperon|s<v[(l + 75) + ^(1 - 7 5 )]&|Ag). 



(4.3) 



Using the interpolating fields [1£ 



Bz(v , s) = u(v, s)(/> v h v , (4.4a) 



B 6 {v, s) = B^v, s)^h v , (4.46) 

where cf> v and ^ are the + and 1 + diquarks, respectively, which combine with the heavy 
quark h v of velocity v to form the appropriate heavy baryon, and the relations [18] 

(O|<V0t|O>= C(vV ), (4.5a) 

(0|#^!|0) = 0, (4.56) 

where ((v ■ v') is a universal baryonic Isgur-Wise function normalized to unity at the zero 
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recoil v ■ v' = 1, we find [] 

<A(t/, ^Iscr^Cl ± 7fe)fe|Ag(i;, s)) = u^l ± l5 )u Ab ((v ■ v'), (4.6a) 

(E°(v', S ')\sa^(l±j 5 )b\A° b (v,s)) = 0. (4.66) 

Therefore, no weak B% — Bq transition can be induced by the b — » sj mechanism in the 
heavy quark limit. So, the first prediction we have is 

T(A° - S° 7 ) « T(A° - A° 7 ). (4.7) 

We will follow Ref. [19] to treat the l/m s corrections to A° — > A°7. First, from the 
relationship between the s quark field and the effective field h^, 



s (x) = e - m ° v - x 



i + — r^- — 



2 2m s 



h { :\ {u 



we get 

iip 



s<V(l ± 75)6 -> J$ 1 1 - ^ J <v(l ± 75)^ 6) . (4.9) 

Applying the result [19] 

(A(t/, sOl^^r^lA?^, s)> = AC(^ • V) V " ~l V v V ' v ] V '" u A (V, s')rTu Ab (v, s), (4.10) 
with the new parameter A being 

A = m A — mb = itl a — m c = m\ — m s « 700 MeV, (4-11) 

we obtain 



(A(v', s')\h$il])a^k v {l ± l5 )h®\A° b (v, s)) 
= hR((v ■ v') u A {v>, s')a, u e^{l ± 75)^ {v, s), 

with h = (— — v ■ v')/(l + v ■ v'), where use has been made of 



(4.12) 



-m K , 



m Ab m Ab 

It follows from Eqs.(4.3), (4.6), (4.9) and (4.12) that 



mA -v' + —k, ft = -U- ( 4 - 13 ) 



a = % ^- 2 F 2 m b V tb V t * s ( 1 + ^ - ^-h) ((v ■ „'), (4.14a) 



In the heavy s quark limit, a A is made of just a strange quark and a scalar diquark with (ud) quantum 
numbers. Contrary to some claims made in the literature, it is not necessary to include a Clebsch-Gordon 
coefficient l/v3 in Eq.(4.6a). 

14 



h =%^ F ™ v » v «( i -^-^M v - v,) - (4 - i46) 

Including QCD corrections gives rise to [20] 

C(v-v') = C(/i)Co(W,/i), (4.15) 

where 

CU) = ( a s(m b ) y 6/25 ( a s {m c ) \-" 27 ( <x a {m 3 ) \ aL ^ 
\a s (m c )J \a s (m s )J \ a s (fi) I 



UJ 



a L {u) = - [^f^j Wo; + vV^l) - lj , (4.166) 

with uj = v ■ v '. There is no obvious choice for the normalization scale in Eq.(4.15). It is 
expected that there will be no large parameters in the function (o(u, /i) if the renormalization 
scale is low [19]. However, since perturbation theory will break down at very low scales, we 
thus choose fi ~ m s so that a s (/i) ~ 1, and C(fi) ~ 1.23 . For the decay A b — > A7, u = 2.63 . 
It is easily seen that the l/m s correction to the A b — ► A7 amplitude is about 50% for 
m s = 510 MeV, which is quite sizeable. This implies that it is important to include higher 
order l/m s corrections. However, this is beyond the scope of the present paper. 

In order to estimate the decay rate of A^ — > A7, we employ two recent models for ((v -v'): 
( (u) = 0.99 exp[- 1.3(a; - 1)] (soliton model [21]), (4.17a) 

/ 2 \ 3 - 5 +i# 
Co(w)= ( J (MIT bag model [22]). (4.176) 

Hence, Co(^ — 2.63) ranges from 0.09 to 0.12 . Substituting (4.14) into (4.2) and using the 
good approximation V tb V t * s = -V cb V* s and the lifetime r(A b ) = 1.07 x 10~ 12 s [15], we find 

B{A° b -► A° 7 ) = 1.31 x 10- 3 |Co(o;)| 2 = (1.2 - 1.9) x 10~ 5 . (4.18) 

Since there is only one strange quark in H& but two strange quarks in Q b and S, it is not clear 
to us how to generalize the above heavy s quark method to the radiative decays S b — > H7 

and Q b — > Q7. 

In the second method, only the 6 quark is treated as heavy. Since 706 = 6 in the static 
limit of the 6 quark, we have the relation 

(A\sia 0i {l ± 75)6|A 6 > = (A|s 7i (l t 7 5 )6|A 6 ). (4.19) 

In the heavy 6-quark limit, there are only two independent form factors in the 7 M and 7^75 
matrix elements [23]: 

(A(p)|s 7 „(l T^)b\A b (v)) = u A (f^ a (v ■ p) +^F 2 AbA (v ■ p)) 7 „(1 TTsK,- (4.20) 
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Recall that F\ — 1 at zero recoil and F 2 = if both A b and A are treated as heavy. The 
relationship between F 12 and the standard form factors defined by 

(A(p)\s 7fl (i - l5 )b\A b ( P ')} = u A ifi(q 2 h» + if2(q 2 K»q u + f 3 (q 2 K 

-(gi(q 2 )j» + m{q 2 )(y^q v + gz{q 2 )q^]u, 

fi{q 2 ) = 9i{q 2 ) = F^q 2 ) 
f2(q 2 ) = 92(q 2 ) = h(q 2 ) 
It follows from (4.3), (4.19) and (4.20) that 



+ ^F 2 (q 2 



gsiq 2 



(4.21) 



(4.22) 



m Ah 



G F e 
72 8^ 
Gf e 
V2 8^ 



F 2 m 6 y t6 v;: 1 



7T 



i^mftV^K 



f.s 



m s 
m b 



,A,A 



A, A/ 



1^(0) - F Ai)A (0) 



A,, A/ 



Ff fcA (0) - F Ai)A (0) 



(4.23a) 
(4.236) 



where F(0) here means -F(g 2 = 0). 



In principle, the Isgur-Wise form factors Fi^{0) can be obtained by first computing the 
form factors /, and g± at zero recoil in a quark model Q and then extrapolating to q 2 = 
under some assumption on their q 2 dependence. In the literature, form factors /, and Oj have 
been calculated for A c —>■ A transition in two different quark models: the nonrelativistic 
quark model and the MIT bag model [24]. The model calculation can be tested by a recent 
CLEO measurement [25] . Assuming a dipole q 2 dependence of the form factors 

*i(0) 



W 



q 2 /ml 



2^2' 



^A C A 



where m* is a pole mass, | the form factor ratio R = F 2 cA / p^ iy i s found to be 



by CLEO for m* 



R = -0.33 ±0.16 ±0.15 

2.11 GeV [25]. From Tables IV and VI of Ref.[24], we find Q 
-0.34 , MIT bag model; 
—0.23 , nonrelativistic quark model. 



R 



(4.24) 



(4.25) 



(4.26) 



5 This is because the quark model calculations are presumably most reliable at zero recoil where both 
baryons are static. 

6 Beyond the heavy quark limit, the pole masses my and ttia for form factors fi and gi respectively are 
different. In general, the pole mass is taken to be the mass of the nearest spin-one meson with the right 
quantum number. For A c — ► A transition, my = ttidji-) = 2.11 GeV, ttia = OId,(i+) = 2.536 GeV, while 
for Ab — > A: my = tubji-) — 5.42 GeV, rriA = mg s (i+) — 5-86 GeV. 

7 Note that the definition of the form factors fi and gi in Ref. [24] follows that in Ref.[26] and is different 
from ours. In terms of our notation, the bag model predictions shown in Table VI for A c — > A form factors 
evaluated at q 2 — with dipole approximation are f\ = 0.46, fi = — 0.19/mA c , g\ — 0.50, gi — — 0.05/mA c , 
while the nonrelativistic quark model predicts f\ = 0.35, fi = — 0.09/mA c , g\ = 0.61, gi = — 0.04/mA c . 
The difference between f h and gi is attributed to l/m c corrections. The form factor ratio calculated in (4.26) 
comes from the form factors f\ and fi- 
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It is clear that while both models' predictions are consistent with experiment, the bag model 
gives a better agreement. 

As for the A& — > A form factors, our strategy is as follows. We will use the bag model to 
compute f± b at zero recoil, and then employ the heavy-flavor-symmetry relation 

Ft\dl c )/Ft%l c ) = F 2 ^ A ( q l b )/Ft bA (ql b ), (4.27) 



where q\ c = m\ c v — q and q\ b = m^ b v — q, together with the CLEO data for R, to determine 
F±2, & t <?m = ( m A 6 — ^a) 2 , which can be extrapolated to q 2 = using Eq.(4.24). In the 
MIT bag model [27], f\ at zero recoil is given by Q 



ft bA (q 2 J = j d 3 r (u b (r)u s (r) + v b (r)v s (r)) , (4.28) 

where u(r), v(r) are respectively the large and small components of the IS1/2 quark spatial 
wave function. Numerically, we find f l b (g^J = 0.95 (see e.g. Ref.[28] for the technique). It 
follows from (4.22), (4.25) and (4.27) that 

Ft b \q 2 J = 1-02, F^\ql) = -0.34, (4.29) 

where only the central values are presented. Since the form factor ratio is measured by 
CLEO using the dipole q 2 dependence, it is natural to employ the same dipole behavior for 
F i b for the sake of consistency. Note that the dipole form 

G(q 2 ) = % ~ ql J™% (4.30) 

(1 - q 2 /mj) 

plays the role of the baryonic Isgur-Wise function £(v ■ v'). Using m* = 5.42 GeV, we find 
G(0) = 0.092, in accordance with £(2.63) = 0.09 ~ 0.12. Putting everything together, we 
finally obtain 

B(A b ^ A + 7) = 0.6 x 10~ 5 . (4.31) 

This model result is close to the previous prediction (4.18). 

We should accentuate the nature of the above two calculations. The first method is 
model independent; but since the effective s quark in the baryon is only of order 500 MeV, 
it is questionable to apply the heavy quark effective theory to hyperons, as evidenced by 
the sizeable l/m s corrections shown above. The second method takes care of the m s effects 
to all orders, but it is hard to assess its reliability. For example, quark model predictions 
for the form factors fi and gi depend on the model content; the form- factor q 2 dependence 
is very sensitive to the choice of the pole behavior: monopole or dipole. Nevertheless, 



i A small and negligible correction to (4.28) is shown in Eq.(5a) of Ref. [24]. 
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we see that both approaches provide a consistent order of magnitude estimate; that is, 
B(A b -> A7) ~ 1 x lO" 5 . 

4.2 Heavy quark symmetry predictions for 2° — > H°7 

Before embarking on quark model calculations for the radiative decays 5° — > S°7, E' c j 
and A° — > £^7, we would like to see what we can learn from applying the heavy quark sym- 
metry to these decays. It turns out that for the antitriplet to antitriplet radiative transition 
S° — » S°7, heavy quark symmetry implies a nontrivial model independent prediction for a/b, 
the ratio of the parity-conserving and parity- violating amplitudes. ^ 

Let us denote 

Oivv = C7 M (1 - 75)6^(1 - 7 5 )w, (4.32a) 

2 »v = C7/,(1 - 7s)« 07,(1 - 75)6, (4.326) 

and apply the interpolating field (4.4a) to the antitriplet heavy baryons to get 

{E° c (v')\Oi^ v \El(v)) = (0\u f (v', s')0„/c„/c^7 jU (l - %)b v d"/ u (l - ^ 5 )ub v (j)lui(v, s)|0) 

= u f (v', s')^f^(l - ^^(v, S )(0|<M 7l ,(l - 75)^110). 

(4.33) 
Lorentz invariance implies that 

<0|^d 7 „(l - j 5 )u(f>l\0) = A(v ■ v')v v + B{v ■ v')v' u . (4.34) 

Therefore, 

(E c (v')\O llxv \E° b (v)) = Ufiv^s'^il-^u^s^Aivv'^ + Bivvy,,}. (4.35) 

Likewise, the matrix element of 2 ^v is 

(E c {v')\O 2fiV \E° b (v)} = u / (^5 / )7M(l-T5)(0|^«#t|0) 7 ,(l-7 5 K^,5). (4.36) 

Again, Lorentz invariance demands that 

(0\(f) v ,ud<j)i\0) = A'(v ■ v')i + B'{v ■ v')f + C'{y ■ v')fj + D(v ■ v'). (4.37) 

Our next task is to recast (4.35) and (4.36) into a more suitable form. Since Oi^ and 
2fll/ are multiplied by F^ v and F^ v [see (2.15)], only the antisymmetric part will contribute. 
Thus we write 

(E° c (v')\O llMU \E° b (v)) = ±u f (v', s')[A{ ltl v v - i„v,) + B( 7 X - 7*<)](1 - !*)<*, *)■ (4.38) 



other decays such as S° — > 3' c 7, A° — > £"7 



We have checked explicitly that heavy quark symmetry alone does not lead to any useful predictions for 

~/0 a 

b ^ " c7j 7V 6 
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By virtue of the equation of motion fu(v, s) = u(v, s), two useful relations can be derived: 

i i 

(l»v u - 7„u M )(l - 75) = -^^(1 + 75) + ^Vt 1 - 7s), (4.39a) 

(7X - 7,<)(1 - 75) = ^v(l - 75) - ^<v(l + 75)/, (4.396) 

which can be further simplified by applying the analog of Eq.(4.13) to the E& — > S c transition. 
As a consequence of (4.33) and (4.39), we obtain 



~b 



{KW)\Oi»A^») = — A U + b m 



«/(*>', «')[<v( 1 + 7s) 



m. 



E 6 



<V(1 -75)]«»(v,«)- 

(4.40) 



As for the matrix element of 02^1/, (4.36) leads to 



(~° c (v')\0 2 , u \~° b (v)) = u f (v', s') [A\ llx f lv - 7^) + B'(j,f lv - lv f%)\ (1 - l5 ) Ui (v, s) 
= iu f (v', s'){A'[(l + f)u^ + (1 - #(V7 5 )] 
+ J BV^[(l + f)-75(l-f)]}^(^s). 

(4.41) 
Applying the analog of (4.13) again to (4.41) leads to 



1—0 



/ , =6 r>' 



(S2(t/)|0^|S»(t;)> = i U' + -^-ff u f (v', s') 



a 



m_ 



jW 



m. 



(l+ 75 ) + _*L(l- 75 ) 



m_ 



Ui{V,S). 

(4.42) 



Finally, substituting (4.40) and (4.42) into (2.15), we obtain 



(E° c (v>)\OF\E°(v)) 



2e 
— 1 T 

1 J 



A + B-f± )t A' + B 



,m— 



x Uf(v' ,s')a ■ F 
+ 



m f 



^c™ &d -r, , ~r \&u m &b m . 



-m c "-m d 



mi „ mi 



' m u rn-b I rn ~ 



(4.43) 






75>tti(U,S). 



From (4.43) we see that although the 3 ^ 3 + 7 transition depends on the unknown param- 
eters A, B, A' and B', a unique tree- level prediction on the ratio of a/b based on heavy 
quark symmetry is nevertheless accomplished: 



c rric rrid \ " m u u inj / m = 

"6 



(4.44) 



This ratio can be tested by measuring the asymmetry parameter a as will be discussed later. 



4.3 Bag model calculations for 2° — > H°(H'°) 7 and A° 



S°7 
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Recall from Sect. II that the effective Hamiltonian responsible for the weak radiative 
decays of heavy baryons is given by 



H cS (bu -> crf 7 ) = ^V cb V: d (c + Ol + C-Oi), (4.45) 

where [cf. Eq.(2.15)] 

0± = d f (F^ + iF^Og + di(F^ - iF^Og, (4.46) 

with Of = Of ± O^, and 

d i =m i [— -)— 2 5", (4.47a) 

d, = ro, (^ - ^) -^-^, (4-476) 

and ?7ii = m& + m u , m/ = m c + m^. Writing 

(^(^OIO^I^^)) = iu f (v',s , )(a lt2 + b 1 ,2'Y5)<T*'u i (v,s), (4.48) 

we get 



(B f (v')\Oi\B t (v)) = -u f (v',s'){d f (a ± + b ± )-d l (a T -b T ) 

+ [d/(a ± + b±) + di(a T - 6 T )] 7 5 }a • F«^, s), 



(4.49) 



where a± = Oi±02, &± = 6i =b 62 , corresponding to the parity-conserving and parity-violating 
matrix elements of Of. It follows from Eqs.(4.45) and (4.1) that 

G F 

a = ^=V cb V* d {c + [df(a + + b+) - d { (a_ - 6_)] +c_[d/(a_ + 6_) - dj(a + - 6 + )]} , (4.50a) 

V2 

O 

6 = ^KbKd {c+[d f (a + + b + ) + di(a_ - &_)] + c_[d/(a_ + 6_) + dj(a + - 6+)]} . (4.506) 

We shall employ the MIT bag model [27] to evaluate the four baryon matrix elements 
a± and b±. Since the quark- model wave functions best resemble the hadronic states in the 
frame where both baryons are static, we thus adopt the static bag approximation for the 
calculation. Because 

u)a xy u\ = 1, u}(T 0z l5 u] =i, (4.51) 

for Vf — Vi — 0, it follows from (4.48) that 

ai, 2 = -i(B f T \0\l\B, t), (4.52a) 

61,2 = ~(Bf t \0%\Bi t). (4.526) 
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Matrix elements 01,2 and bi^ in the MIT bag model can be expressed in terms of four- 
quark overlap bag integrals (see e.g. Ref. [28] for the method). To do this we need the 
spin-color wave functions of the baryons involved such as 

~° = j={(bus - bsu) XA + (12) + (13)}, 
S o = _^[ c ^ Xs + (i2) + (13)], 

where abcx s — (2a^c^ — a)b^c^ — a^b^rf)/\/6, abcxA — (cfitfc^ — <rfb^rf)/\/2, and (ij) means 

permutation for the quark in place % with the quark in place j. As an example of the bag 

model evaluation, we look at the decay S° — > S°7. With the wave functions given by (4.53), 

we find 

(E° c T \blhbblb2ual\E , T) = \, (4.54a) 

6 

(S° T \bl c b u bi d b2u<7 z 2 \Z b T) = 0, (4.546) 

(S° T \b\ c b lb bl d b2u(cri + <J2- ~ (7i_<7 2+ ) |E° T) = 0, (4.54c) 

(£° T IblhubiM^ T) = ^ (4.54d) 

(S° T \blb lu bi d b 2b a z 2 \~° b T) = i, (4.54e) 

(£° T |&LMLM*i+*2- " ^i-^2+) |2g T) = ^- (4.54/) 

For the decay S° — > 2^7, we find finally 

4tt r-R^.2 



oi = -6 2 = -f J r 2 dr(v c u b + u c v b )(v d u u - u d v u ), 



b\ = -a 2 = — f 1 J r 2 dr[(3« c w fe - v c v b )(u d u u + v d v u ) + (w c v b - v c u b )(u d v u - i> d u u )]. 

(4.55) 
Consequently, 

a± = ^±, 6± = T^±, (4-56) 

with 

47T r R 

/+=—-/ r 2 dr[(u c v b + 3v c u b )(v d u u - u d v u ) + (3w c w fe - v c v b )(u u u d + i> u i> d )], (4.57a) 
3 Jo 

4tt f R 
I- = — r 2 dr[(3u c v b + v c u b )(v d u u - u d v u ) - (3u c u b - v c v b )(u u u d + v u v d )\. (4.576) 
3 Jo 

From Eqs.(4.50) and (4.56) we obtain 

a = -^V cb V: d c4d f I_ - dil+), (4.58a) 

b = -^V cb V: d c.(d f I. + dj + ). (4.586) 
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Several remarks are in order, (i) We have explicitly confirmed that the operator 0+ does 
not contribute to the baryon transition matrix elements as the baryon-color wave function is 
totally antisymmetric. This is ascribed to the fact that 0+ is symmetric in color indices, as 
can be seen by applying the Fierz transformation to the effective operator for the amplitude 
(2.10) and by noting that the photon interaction is color singlet. From Eq.(4.43) we conclude 
that A = 4A' and B = AB'. (ii) In the isospin limit we have I + = — J_, so there is only one 
independent bag integral for the decay amplitude of S b — > S°7. In the same limit, we find 
that 



mt m 

e c — — e d 



d \ u m u "rrtbj 



^b —c I 



rfl f Tn f ( m» m» 

c m c " rrid \ u m u " mi, 



(4.59) 



Comparing this with (4.44), it appears that the bag model does not predict correctly the 
ratio a/b. This seeming inconsistency comes from the static bag approximation we have 
adopted. In the rest frame of the initial baryon, one can show that the ratio of the masses 
mi and rrif is given by 

rrii 



r = 



1 + Vf/c 

m f \l-v f /c 



(4.60) 



The ratio a/b is in principle a function of the quark masses, the bag parameters and the 
velocity Vf. But in the static bag approximation, we always have r — 1. In order to get the 
heavy quark symmetry prediction for a/b, we thus need to utilize a moving bag to describe 
the recoil effect of the final baryon state. The net effect should be that (4.58) is modified to 
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(4.61a) 



(4.616) 



as implied by Eq.(4.43). Later we will use (4.61) rather than (4.58) to compute the decay 
rate and branching ratio for the decay S° — ► S°7. (iii) Experimentally, the ratio a/b can be 
determined by measuring the asymmetry parameter a when the initial baryon is polarized 
with the polarization vector sf 



dT(B t -> B n ) 

dn 



r/ iB - 



Bfy)(l + asi-p f ), 



where 



a 



2Re(a*o) 

|a| 2 + |6| 2 ' 



(4.62) 
(4.63) 



For completeness, we shall write down the results for the remaining two decay modes of 
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the bottom baryon. For S° — > E' c 7, we get 

0l = -6 2 = ^ Jjf r 2 dr[u c u b (3u u u d - w u w d ) + 2v c u b (u u v d + w u u d ) - v c v b (u u u d + <vu d )], 
61 = -02 = -^ Jcf r 2 dr[v c u b {hu u v d - v u u d ) - u c v b {3v u u d + w M v d ) 

+ [u c U b + V c V fc )(3w„W rf - V u V d )\. 

(4.64) 



As a result, 

with 

4vr /''< 



a± = 7^4, &± = T~4, (4.65) 



1+ =-- -rrjn I r 2 dr[3(u c u d + v c v d )(3u u u b - v u v b ) + (3u b v u + v b u u )(u d v c - v d u c )}, (4.66a) 

Air f R 
I'_ = —y= / r 2 dr[{3u c u d - v c v d ){u u u b + v u v b ) + (u b v u - v b u u )(5u d v c - v d u c )}. (4.666) 

The resulting amplitudes for the decay S° — > S' c 7 are 

a = -^VaKfi-idflL - dj' + ), (4.67a) 

As for the transition A° — > £°7, we find that its amplitude is the same as that of H° — > S' c 7 
except for a different overall normalization factor. More precisely, 

o(Ag - S° T ) = V2a(H° - H'° 7 ), (4.68a) 

&(Ag-E° 7 ) = V^g-S^). (4-686) 

Finally, we come to numerical estimates. For the bag parameters we use [27] F] 

m M = m d = 0, m s = 0.279 GeV, m c = 1.551 GeV, m fe = 5.0 GeV, 
x u = 2.043, a;, = 2.488, x c = 2.948 , x b = 3.079 , i* = 5.0 GeV" 1 , 

where i? is the bag radius. The eigenvalues x q 's are determined by the transcendental 

equation 

x 

tanx= i- m R-( X 2 + m m 2 y/f (4 ' 70) 

Numerically, the relevant bag integrals are found to be 

1+ = -I- = 2-443 x 10- 3 GeV 3 , 

I' = 3.720 x 10~ 3 GeV 3 , I' = 1.267 x 10~ 3 GeV 3 . '' 



10 It should be stressed that except for the bag quark masses used in (4.69), all the light quark masses 
employed in the present paper are of the constituent type [see Eq.(3.10)], as explained in Sect. II. 
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Putting everything together and using m s = 5809 MeV [29], m , = 2573 MeV, we finally 
obtain the decay rates 

T(E° b -► H° 7 ) = 3.95 x 10~ 20 GeV, (4.72a) 

T(S° -► E'° 7 ) = 3.54 x KT 19 GeV, (4.726) 

T(A° b -► S° 7 ) = 6.65 x KT 19 GeV, (4.72c) 

the branching ratios 

£(~° - ~° 7 ) = 6.4 x 10- 8 , (4.73a) 

B(~° b -> E'° 7 ) = 5.7 x KT 7 , (4.736) 

£(A° -► S° 7 ) = 1.1 x 10- 6 , (4.73c) 

for r(E°) ~ t(A°) = 1.07 x 10~ 12 s [15], and the decay asymmetry 

a(~°^~° 7 ) = -0.47, (4.74a) 

a(S°^~'° 7 ) = -0.98, (4.746) 

a(Ag->E°7) = -0.98. (4.74c) 

Note that the prediction of a for the decay mode 5° — * S°7 is based on heavy quark 
symmetry [cf. Eq.(4.44)]. Since the baryonic matrix elements are evaluated under the static 
bag approximation which amounts to a maximal overlap of wave functions, the decay rates 
and branching ratios given by (4.72) and (4.73) for 2° — ► S' c 7 and A° — * S°7 ought to be 
regarded as the most optimistic estimates, and the respective decay asymmetry parameters as 
order of magnitude estimate. Nevertheless, the sign of a given in (4.74) is more trustworthy. 
Finally, a comparsion of (4.73) with (4.18) leads to the conclusion that the weak radiative 
decays of bottom baryons are indeed dominated by the electromagnetic penguin mechanism. 

V. Conclusions 

Nonpenguin weak radiative decays of heavy hadrons are characterized by emission of a 
hard photon and the presence of a highly virtual intermediate quark between the electro- 
magnetic and weak vertices. We have argued that these features should make possible to 
analyze these processes by perturbative QCD. 

In this work we have found in tree approximation that these processes are describable by 
an effective local and gauge invariant Lagrangian. This Lagrangian leads to a unique heavy 
quark symmetry prediction for the asymmetry parameter in the decay S° — * S°7. Other 
interesting results are obtained by making use of the effective Lagrangian in conjunction with 
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the factorization approximation for heavy meson decays and the MIT bag model for heavy 
baryon decays. In particular, the branching ratio for B° — > D *'-/ is found to be 0.9 x 10 -6 . 
This is very important for the experimental interpretation of the inclusive measurement of 
B° — > 7+ anything and its relation to the penguin dominated decay B° — > K°* / -f. In Section 
4 we have presented two different methods for estimating the rates of A& — ► A7. Both 
approaches provide a consistent order of magnitude estimate, namely B(Ab — > A7) ~ 1 x 10~ 5 . 
We conclude that weak radiative decays of bottom hadrons are dominated by the short- 
distance b — > s'-f mechanism. 

The factorization method has been known to be reliable for nonleptonic decays of heavy 
mesons. So the prediction for B° — > D* 0/ ~f based on this method should also be reliable. 
For the heavy baryon sector, we have to resort to the static bag approximation even though 
the initial and final heavy baryons move with substantially different velocities. This is a 
serious drawback in our work. The bag results obtained in Section 4 can only be regarded 
as order of magnitude's rough estimates. It remains an important theoretical question how 
to incorporate the relative motion between two bags in the MIT bag model. 

Cabibbo-allowed weak radiative decays of charmed hadrons give rise to strange hadron 
in the final states. The constitutent s quark, whose mass is of order 500 MeV, is not very 
heavy. We can make rough estimates of the branching ratios for the aforementioned decays 
by treating the s quark as heavy and applying the formalism developed in this paper. Even 
though the l/m s corrections are expected to be (50 — 100)%, they should not alter the order 
of magnitudes. In this way, we find 

B{D° -> iT°7) = 1.1 x 10" 4 , (5.1) 

£(A+ -> S+ 7 ) = 4.9 x 10~ 5 , i3(S°^S° 7 ) = 3.1 x 10" 5 , (5.2) 

a(A+ -► S+ 7 ) = -0.86 , a(E° c -► S° 7 ) = -0.86 , (5.3) 

where we have used the lifetimes for D, A c , S c from PDG [15]. The branching ratio for 
D° — > K*° + 7 is quite sizeable. These decays should be searched for experimentally. 

Finally, we observe that in the weak radiative decays of bottom hadrons the highly 
virtual quark's squared invariant mass is of order ml or smaller. It is therefore appropriate 
to employ the renormalization group improved weak interaction Hamiltonian (2.4) with a 
renormalization scale \x = nib as a starting point. If one wishes to use a renormalization 
scale smaller than mb, one must reanalyze the one loop corrections in a heavy quark effective 
theory including diagrams with the photon inside the loop. However, we have not considered 
such an analysis. 
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Figure Captions 

Fig. 1. IF-exchange diagrams contributing to the quark-quark bremsstrahlung process 
b + d — >c + m + 7 induced by the four-quark operator Oa defined in Eq.(2.6). 

Fig. 2. Same as Fig. 1 except for the operator Ob defined in Eq.(2.6). 

Fig. 3. Possible pole diagrams contributing to the radiative decay B° — > D* 0/ ~f. 
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